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' Selected applications of symmetry methods in the atmospheric sciences are 

I . reviewed briefly. In particular, focus is put on the utilisation of the classical 

Lie symmetry approach to derive classes of exact solutions from atmospheric 
models. This is illustrated with the barotropic vorticity equation. Moreover, 
the possibility for construction of partially-invariant solutions is discussed for 
this model. A further point is a discussion of using symmetries for relating 
different classes of differential equations. This is illustrated with the spherical 
and the potential vorticity equation. Finally, discrete symmetries are used to 
1-^^ ' derive the minimal finite-mode version of the vorticity equation first discussed 

, by E. Lorenz (1960) in a sound mathematical fashion. 
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>. ■ 1 Introduction 

(N 

Dynamic meteorology is concerned with the mathematical theory of atmospheric 
T^lj- ■ motion. It lies the corner stone for daily weather prediction as it prepares the 

, grounds for numerical computer models, without which reliable forecasts are 

I hardly imaginable. However, since the advent of capable supercomputers and 

the accompanying shift towards methods for exhausting their capacities, analyti- 
cal investigations of the governing equations have somewhat taken a back seat. 
However, as evident, numerical models need benchmark tests to check for their 
^ ' reliability and hence exact solutions of the underlying mathematical models are 

■ still of great value. Testing whether a forecast model is able to reproduce a known 

exact solution of the original equation may serve as a first consistency check. It 
follows that techniques for obtaining such solutions in a systematic way are of 
rather pretty importance. For this purpose, the classical Lie symmetry methods 
are well-suited. 

Lots of dynamical models in use in the atmospheric sciences are adapted forms 
of the Navier-Stokes or the ideal Euler equations, taking into account both the 
rotation of the earth and the anisotropy of the atmosphere (the region of interest 
for weather prediction extends about 10 km in the vertical but several thousands 
of kilometers in the horizontal direction). For large-scale dynamics (i.e., horizontal 
length scale about some 1000 km) or for sake of conceptual simplicity, it is possible 
to restrict oneself to two-dimensional models. The most relevant example of such 
a model in dynamic meteorology is the barotropic vorticity equation. It is derived 
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from the incompressible Euler equations in a rotating reference frame by using 
the stream function (or vorticity) as dynamic variable. 

In this contribution, we discuss symmetries and group-invariant solutions of the 
barotropic vorticity equation. Moreover, the construction of partially-invariant so- 
lutions is shown. These issues are addressed in section [21 Subsequently, in section 
[3]the usage of symmetries for finding related differential equations is demonstrated 
using the potential and the spherical vorticity equation. It is shown that in both 
cases the rotational term in the equation can be canceled using suitable point 
transformations. In section [4] discrete symmetries of the vorticity equation are 
used to systematically re-derive the Lorenz (1960) model. This contribution ends 
with a short summary and discussion of future research plans, which can be found 
in the final section \5\ 

2 Group-invariant solutions and barotropic vorticity 
equation 

It is worth considering the Lie symmetry problem of the barotropic vorticity equa- 
tion since to the best of our knowledge, this equation has not been investigated 
thoroughly in the light of symmetries before. Either only the symmetries and 
some exact solutions were computed without reference to classification of inequiv- 
alent subgroups [71 18] or the classification itself was not done in most complete 
fashion |JJ. 

Using the stream function-vorticity notation, the barotropic vorticity equation 
in Cartesian coordinates reads: 

Ct + IpxCy - IpyCx + (3lpx = 0, C = Ipxx+'Ipyy, (1) 

where C stands for the vorticity, ip is the stream function and /3 = const is a 
parameter controlling the north-south variation of the earth angular rotation. It 
can be expressed as the north-south change of the vertical Coriolis parameter 
/ = 2Qsm(p via /3 = df/dy, where is the absolute value of the earth angular 
rotation vector and if denotes the latitude. 

Equation ^ admits the infinite-dimensional maximal Lie invariance alge- 
bra generated by the operators 

V = tdt - xdx - ydy - Sipd^, dt, dy, 
X{f) = fmx - f{t)yd^, Z{g) = g{t)d^, 

where /i, / and g are arbitrary real-valued time-dependent functions. For a sys- 
tematic group-invariant reduction of the vorticity equation by means of using sub- 
algebras of B"^ , it is necessary to compute the corresponding optimal system of 
inequivalent subalgebras first [lOj. For one-dimensional subalgebras it reads [HH] 



(P), {dt + cdy), {dy+X{f)), {X{f)+Z{g)), 
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where c = {0, ±1}. For two-dimensional subalgebras the optimal system is [1] 

{V,dt), {V, dy + aX{l)), {V, X{\tn+cZ{\tr^)), {V,Z{\tr^)), 
{dt + bdy, X{e^^) + Z{{aht + c)e'^*)> , (5* + hdy, Z{{abt + c)e''')) , 
{dy + X{f'), + Z{g^)) , (a, + Xif), Z{g^)) , 

{X{f')+Z{g'), X{f)+Z{g^)), 

with a,b,c being arbitrary constants and /* and g^, i = 1,2, being arbitrary 
real-valued functions of time. In the fifth and sixth subalgebra, additionally the 
condition abc = has to hold. In the last subalgebra the pair of functions {f^,g^) 
and has to be linear independent. Using the above two-dimensional sub- 

algebras allows to reduce ([T]) to ODEs. Moreover, it is straightforward to show 
that the third and the fourth case of one-dimensional subalgebras lead to com- 
pletely integrable PDEs. Hence, since all but the first two-dimensional subalgebra 
are extensions of special forms of the third or fourth one-dimensional subalgebras. 
That is, only reduction by means of the first two-dimensional subalgebra may give 
an essential new result. 

Based on the above classification, it is possible to derive classes of inequivalent 
group-invariant solutions fl]. As a notable example, we investigate reduction by 
means of {dy + X{f)). The invariant functions of this subalgebra are 

1 , 2 

P = x-fy, q = t, v = 'ilj + -fy, 

which allows one to reduce ([I]) to the classical Klein-Gordon equation Vpq+fiv = 0, 
where 

dg . . f" ^h{q) ({l + f)f"y 
p = p, V = V —p H — h 
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For the meteorological application, the case / = const and making a harmonic 
ansatz for v is the most relevant: The corresponding solution represents the clas- 
sical Rossby wave, which to a large extend governs the weather regimes in the 
mid-latitudes. 

Considering ([1]) as a system of two equations using Tp and C as dependent vari- 
ables, it is possible to compute partially-invariant solutions [TI]. This should be 
illustrated using the subalgebra Z{g)). Due to the presence of the basis 

element Z(g), no ansatz for ip can be chosen. However, it is still possible to make 
an ansatz for C. Because of the joint invariance under A'(l) and Z(g) we have 
C = C(^) y) and thus ([l]) transforms to 

Ct + TpxiCy + P) = 0, C = '>Pxx+^yy 

Introducing the absolute vorticity tj = ( + (3y we have to distinguish two cases. 
For r]y = the solution of the vorticiy equation is 

ilj = ^{t,x,y) -]:Py^ + lr]{t)y^, 
b 2 
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where satisfies the Laplace equation "^^x + ^ yy = 0. For i]y ^ the solution 
reads 

where lo = g^y + and g^,g^, f^, are all real- valued functions of time. 

3 Symmetries and the ineffective earth rotation 

3.1 The spherical vorticity equation 

Although the vorticity equation in Cartesian coordinates allows one to study some 
prominent features of large-scale geophysical fluid dynamics, it does not take into 
account for the earth's sphericity. For this purpose, it is necessary to study 
vorticity dynamics in a rotating spherical coordinate system. The corresponding 
equation then is [13]: 

1 20, 

Ct + — (V'aCm - V'mCa) + —i^x = 0, (2) 

where ■0 is the (spherical) stream function and C the (spherical) vorticity. 



c = 4 



^-^Vaa + ((1-/^')V'm), 



Rather than using A (longitude) and ip (latitude) as spatial variables, it is advan- 
tageous to use A and fi = sin 99. The mean radius of the earth is denoted by a. 
The maximal Lie invariance algebra of ([2]) is generated by the basis operators 

V = tdt - {ip - On)d^ - fttdx, dt, Z{g) = g{t)d^, Ji = dx, 

J2 = ^^ n ^ ^A + /- ^ ((1 - f^ )d^ + ^9^) , 

J3 = \ .^ dx V ((1 - i^^)Q^ + ^d^) , 

where g runs through the set of smooth functions of t. It is straightforward to 
map the algebra gf^ with 7^ to the algebra go by means of the transformation 

t = t, /i = /i, A = A + ^It, Tp = Tp — Qfi. 

Moreover, this transformation also maps the equation ^ with 7^ to the 
equation of the same form with f2 = 0. That is, it is possible to disregard the 
rotational term by setting 17 = for all practical calculations and finally obtain 
the corresponding results for the case $7 7^ by applying the above transformation. 
We note in passing that this transformation was already used by Platzman to 
transform the vorticity equation to a zero angular momentum coordinate system. 

For optimal systems of one- and two-dimensional subalgebras of go and the 
computation of group-invariant solutions of ([2|) in a fashion similar to the previous 
section, see 1:4 1. 
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3.2 The potential vorticity equation 

An extension of the classical barotropic vorticity equation is given through the 
potential vorticity equation. For the flat topography, the equation reads 

Ct - Flpt + ipxCy - i'yCx + f^tpx = 0, ^ = + iPyy, (3) 

where F is the ratio of the characteristic length scale to the Rossby radius of 
deformation [12]. In the Lagrangian view, the above equation can be understood 
as individual conservation of the potential vorticity q = C + f — Fip. In this 
light, the extension to the barotropic vorticity equation, which in turn states the 
individual conservation of absolute vorticity r] = C + /) becomes most obvious: 
While barotropic dynamics takes place solely in two dimensions, the additional 
term —Fip accounts for a variation of the fluid height in the vertical. In this 
respect, the potential vorticity equation is especially suited for shallow water 
theory. 

Determining the symmetries for the case /? = and P ^ shows that both 
algebras and hence also the two corresponding equations are mapped to each other 
by applying the transformation [3] 

t = t, x = x + —t, y = y, ip = ip-—y. 

Again, this shows that there are models in the atmospheric sciences, in which the 
rotational terms are apparently not of prime importance. 



4 Symmetries and finite-mode models 

There is a long history in dynamic meteorology to convert the governing nonlinear 
PDEs to systems of coupled ODEs by means of series expansions, followed by a 
reasonable truncation of the series in use. Among one of the first models that 
was analyzed in this way again was the barotropic vorticity equation in Cartesian 
coordinates. This was done by Lorenz [9j in an ad-hoc fashion: He first of all 
expanded the vorticity in a a double Fourier series on the torus 

C = ^Cmexp(im • x), 

m 

where x = xi + yj, m = mii + m2j, m = mik\ + 771-2^, k = const, / = const, 
i = (1, 0, 0)"'", j = (0, 1, 0)^, mi and m2 run through the integers and the coefficient 
Coo vanishes. Afterwards he substituted this expansion in ([T|) for /3 = yielding 
the spectral form of the vorticity equation: 

^ = - L (k • [m X m]) , (4) 

where k = (0,0, l)""". Finally, he restricted the indices of the Fourier coefficients 
Cm = l/2(^m — iB^) to only run through the indices { — 1,0, 1}, leading to an 
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eight-component initial model. He then noted that if the imaginary parts i?m of 
the coefficients at hand vanish initially, they will remain zero for all times. More- 
over, if the real coefficient An = — _i initially, this equality is also preserved 
under the finite- mode dynamics. These observations allow to reduce (jl]) to the 
following three-component model 



Now it is instructive to see, whether it is possible to give a more rigorous 
justification of the two observations by Lorenz [2]. The key for this investigation 
is to start with the generators of admitted mirror symmetries of ([T]) for /? = 

ei: {x,y,t,'il;) ^ {x,-y,t,-^l)), 
e-2- {x,y,t,i)) ^ {-x,y,t,-'il)), 
63 : (x, y, t, -0) ^ (x, y, -t, -Tp). 

and induce them in the space of Fourier coefficients by means of series expansion. 
The induced symmetries are 



ei : 


Cmim2 


Cmi,— 1712 7 


62: 


Cmim2 




63: 


Cmim2 


Cmim2 ; 



t ^ -t. 



Moreover, the translations with the values ir/k and n/l in directions of x and y, 
respectively, induce proper transformations of the Fourier coefficients: 



The transformations ei, 62, p and q act on dependent variable and thus may be 
useful for reducing the number of Fourier coefficients at hand. Therefore, the 
finite dimensional symmetry group that is relevant for us is generated by these 
four elements and has the structure G ~ Z2 © Z2 © © Z2. Selecting various 
subgroups of G allows to reduce the eight-component initial model to fife-, four- 
and three-component submodels, respectively. The three-component submodel 
by Lorenz is derived upon using the subgroup S = {I,pqei,pqe2,eie2}- The 
transformation 6162 accounts for his first observation, i.e. Bmim2 = 0, while the 
transformation pqei yields the identification An = —Ai-i, which justifies the 
second observation. In this way ([5]) is derived merely using symmetry techniques. 

This example shows that it is possible to derive a consistent finite-mode model 
by means of a sound mathematical method. Having such a method at ones disposal 
is especially useful since up to now there are only few criteria for the selection of 
modes in finite-mode models. 




(5) 



where, A 



Aoi.F := Aia.G := Ai^_i. 
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5 Conclusion and future plans 

In this paper, we have reviewed some possible apphcations of weU-estabHshed 
symmetry methods in the atmospheric sciences. The technique of constructing 
group-invariant solutions allows to systematically re-derive the well-known Rossby 
wave solution of the barotropic vorticity equation. Moreover, relating differential 
equations by mapping their Lie algebras to each other enables one to cancel the 
terms due to the earth's angular rotation in both the potential and spherical vor- 
ticity equation. Form the physical viewpoint, this result is somewhat astonishing 
since these rotational terms are in fact to a great extend responsible for the overall 
complexity in formation of weather pattern. 

In future works, the classical method of Lie reduction should be applied to 
more sophisticated models of geophysical fluid dynamics, including the quasi- 
geostrophic model [6| and some classical convection model |5|. Moreover, the usage 
of symmetries for a systematic determination of conservation laws in atmospheric 
science should be examined. This again has potential application in providing 
consistency checks for numerical integration schemes. 
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